Quantum Dynamics of Multiferroic Helimagnets: a Schwinger-Boson Approach 
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We study the quantum dynamics/fluctuation of the cycloidal helical magnet in terms of the 
Schwinger boson approach. In sharp contrast to the classical fluctuation, the quantum fluctuation 
is collinear in nature which gives rise to the collinear spin density wave state slightly above the 
helical cycloidal state as the temperature is lowered. Physical properties such as the reduced elliptic 
ratio of the spiral, the neutron scattering and infrared absorption spectra are discussed from this 
viewpoint with the possible relevance to the quasi-one dimensional LiCu202 and LiCuVCU. 
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Frustration, competition between interactions, in mag- 
nets has been an intriguing issue in the field of clas- 
sical/quantum magnetism over several decades. In the 
usual case, even with the competing exchange interac- 
tions Jjj's, their Fourier transform J(q) has the maxi- 
mum at some wavevector q = Q, and the classical ground 
state becomes the helimagnet This is because of the 
constraint of the fixed length of the classical spin, i.e., 
\Sj\ =fixcd. In strongly frustrated quantum magnets, on 
the other hand, the long-range order is possibly destroyed 
and novel ground states without magnetic order may be 
realized. Many possibilities such as chiral spin liquid [2L 
spin-nematic [3j and spin-Peier Is/ valence-bond-crystal [4] 
states are theoretically proposed. Another possibility is 
a magnetically ordered state realized by the order-by- 
disorder mechanism when the corresponding classical sys- 
tem has continuously degenerate ground states [Bj]. 

Recently a renewed interest has been focused on the cy- 
cloidal helimagnets from the viewpoint of multiferroics, 
which exhibit both magnetic and ferroelectric properties 
[U 0] . These materials shed some new light on the frus- 
trated magnets since the electric polarization is cl osely 
related to the vector spin chirality Si x Sj [j.lflfllillllfia]. 
which has been the subject of intensive interests. Namely, 
it was found that the electric polarization(P) produced 
by the neighboring spins (Si and Sj) can be written as 



P = aeij x (Si x Sj), 



(1) 



where denotes the unit vector connecting the sites 
i and j. This relation has a physical interpretation in 
terms of spin current induced between noncollinear spins 
due to frustration [8|. 

Magnetic materials with the finite vector spin chi- 
rality include wide range of systems such as three di- 
mensional(3D) magnets i?Mn0 3 (R =Gd, Tb, Dy) with 
spin S = 2 [3, 0, EH, Ell, the kagome staircase com- 
pound M3V2O8 with 5* = 1 [l7j |. 5=1/2 quantum spin 



chains LiCu 2 2 [H, [3, LiCuV0 4 US and the quasi- 
one-dimensional(lD) molecular helimagnet with S = 7/2 
2l| . Depending on the temperature, dimensionality, 
and the magnitude of the spin S, the role of the classi- 
cal/quantum spin fluctuations differs and the theoretical 
studies on these fluctuations are needed for the consis- 
tent interpretation of the phase diagram and the phys- 
ical properties of these systems. Especially, the low di- 
mensionality enhances both thermal and quantum fluc- 
tuations leading to the breakdown of the conventional 
(classical + spin wave) picture for helimagnets. The 
possible chiral spin states without the magnetic long 
range ordering have been proposed theoretically for clas- 
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241 ] and quantum [25|, |26|, |27[ spin systems. 
However, the systematic study of the quantum fluctua- 
tion in the helimagnets including the finite temperature 
effect is rare, which is addressed in this paper and will 
be complementary to the works mentioned above. 

In this paper, we study the quantum/thermal fluctua- 
tion in the helimagnet in terms of the Schwinger Boson 
(SB) approach. The advantage of the SB method is that 
it can describe the length of the ordered moment as a soft 
variable. Namely, in the constraint on the Schwinger bo- 
son number at each site, 



(2) 



it can be decomposed into the condensed (classical) part 
and the fluctuating part. Therefore, the degrees of clas- 
sical/quantum fluctuation and the ordered moment can 
be described in a unified fashion in this method [28|. In 
the SB language, the paramagnetic to collinear transi- 
tion is described by the density wave instability of bosons, 
while the collinear to helical one corresponds to the Bose- 
Einstein condensation(BEC) of SB. 

Effective model — We study quasi- ID and two- 
dimensional (2D) Heisenberg models with the exchange 
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FIG. 1: Schematic lattice structure and exchange interactions 
of the effective spin model. Ji is a ferromagnetic, while J2 is 
antiferromagnetic interactions. 11/2-coordinates and afcc-axes 
are also shown. 



interactions shown in Fig. [T] where J\ is ferromagnetic, 
while J 2 are antiferromagnetic, leading to the frustration. 
The interchain/interplane interaction J± is assumed to 
be sufficiently weak, and will be treated by the mean 
field theory. The spin-S 1 operators are represented by 
SB as S a = b\{o^ a , /2)b a > , where a a (a = x,y,z) are 
the Pauli matrices and the repeated indices are summed 
over. First, we assume that the resonating- valence- 
bond(RVB) correlation is dominant and neglect the other 
mean-field decoupling. This assumption is valid for the 
low-dimensional multiferroic s l29ll such as LiCuV04 [2(| , 
LiCu 2 2 [3] and Ni 3 V 2 8 [17]. The mean-field Hamil- 
tonians of the quasi-lD model is given by 

Hf 1 ! = X)(A - 277! cos k x )bl a b ka 

ka 

+ ^2 2 I ? 72 sin^fe^) + rj±(smk y + sin k z )]b\^_ ki + h.c. 



+ 2M(r 1 2 1 /J 1 + V l/J2 + 2r 1 l/J ± -S\), 



(3) 



where M is the total number of sites and b^ a is the Fourier 
transform defined by 6 Jcr = J2 k e~ lk ' Rj bk a /VJ7. In 
, A denotes the chemical potential for the bosons and 
the order parameters 771, 77 2 and rj± are Ji{bj (r bi+£ t<T )/2, 
J2{bi V e v<J b i+ -2x,a) I '(2i) and J±{bi V e v „b i+ ^ (T )/{2i) (e = 
y, z), respectively, with ej.f = — e-fj. = 1- RVB order 
parameters are assumed to be real and spatially uni- 
form. In a parallel way, we can derive the quasi-2D 
mean- field Hamiltonian ■ The Hamiltonian 
can be diagonalized by the Bogoliubov transformation as 
#i M d F = E feCT ^(fe)(7L7w+l/2)-2^A(S'+l/2)+con S t, 
with the dispersion relation oj(k) 2 = (A — 2r/i cosk x ) 2 — 
(2?7 2 sin(2/c :I ;)-|-277^(sinfcy + sinfc 2 )) 2 . The transformation 
between 7fc CT and bk a is given by 



bk-\ 



-hi 



cosh 9k sinh 9k 
sinh 9k cosh 9k 



Ik] 
1-hi 



(4) 



withtanh20fc = — (2r7 2 smk x + 2ryj_ (sin fc y +sinA: z ))/(A — 
2771 cos k x ). The chemical potential A is determined by 
the condition ([2|). 77's are obtained by minimizing the 
mean- field free energy F MF . Figure [2] shows the numeri- 
cally obtained 771, ?7 2 and the gap A(T) = u>(Q/2) of the 
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FIG. 2: RVB order parameters 771 and 772 and the gap A of 
the S =1/2 ID model with varying J1/J2 at zero tempera- 
ture. Inset shows the temperature dependence of 771 and 772 
at J1/J2 — 0.5. We use the unit J2 — fee = 1. 
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FIG. 3: a) Schematic energy dispersion of 7-particles. Minima 
are at k = ±Q/2. b) The reorganization of the SB due to 
the collinear magnetic order. The origin of the momentum 
k is shifted by ±Q/2. c) Goldstone and amplitude modes 
associated with the BEC in the helical phase. 



ID spin-1/2 model as a function of Ji/J 2 [3(j. The tran- 
sition temperature of ?7 2 is analytically given by Trvb = 
J 2 (S+ 1 /2) / ln(l + 1/S). We have also numerically stud- 
ied the 2D model at finite temperature and obtained 
similar results. From 77's, we can estimate the minima 
of the dispersion w(fe) as ±Q/2 = ±(Q/2, tt/2, tt/2). Q 
is determined to satisfy (A — 2?7i cos(<3/2))?7i sin(Q/2) = 
4(772 sin Q + 2tj±)t]2 cos Q. 

To describe the low-energy physics of the model, it is 
useful to construct an effective continuum model. First, 
we suppose that 77's are non-zero. Next we expand the 
dispersion around the minima up to quadratic order in 
fe± Q/2. The effective dispersion relations of 7-particles 
are those of massive relativistic bosons and explicitly 

given by f2(fe) 



2 , where fen 



/A(T) 2 + C 2 |fc|||" -f Cj_\kx\~i wnere |S|| 

is the vector along (within) the chain (plane) while k± is 
that perpendicular to the chain (plane). The spin wave 
velocities C|| and c± can be written in terms of 77's, in 
principle. Now the effective Hamiltonian of our system 
is 



Slka 



1/2), 



(5) 



ka a=± 



where a = + (— ) indicates that the momentum is around 
+Q/2 {-Q/2), When the gap A(T) = u(Q/2) vanishes, 
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f2's are the linear dispersions of the Goldstone modes. 

Collinear phase — To study the instability toward the 
magnetic ordering, we consider the mean-field decoupling 
of the interchain/interplane interaction corresponding to 
the density wave formation of the SB and treat the re- 
sulting one/ two-dimensional problem [3l|, (32|. The total 

H int with 



_ rrMF 
— n \Dj2D 



hamiltonian is given by H lD j 2D 

H int = zJ x {\a + ib\ 2 - [(a - ib) ■ S Q + h.c.]}, (6) 

where z is the coordination number along 
interchain/interplane direction and Sq = 
Sfe b\. + Q a (cr aa i /2)b ka > . Here a and b are mean 
fields for (Sq) = a + ib, and collinear and helical orders 
are expressed by them as (a x b) = and (a x b) ^ 0, 
respectively [24|. The interaction between 7-bosons, 
when translated from that between 6-bosons by Eq.(|3}, 
is enhanced near the bottom of the dispersion, inversely 
proportional to the gap A(T) in Fig. 3. a, inevitably 
leading to the density wave instability before the 
occurence of BEC. From the rotational symmetry in 
spin space, we can set a z = b z = without loss of 
generality. By introducing s = (a x — ia v ) + i(b x — ib y ) 
and t = (a x — ia v ) — i(b x — ib v ), we can rewrite H mt 

as H** ~ (zJj2)Z^ {\s\ 2 ~ K-Q/ 2T W2i + 
h.c.)} + (zJ ± /2) E^oiltl 2 ~ K+Q/2A-Q/24 + h - c -)}- 
The summations over k are restricted to around since 
our continuum model is valid only in the low-energy 
region. The free-energy density corresponding to the 
Hamiltonian H = H^^ 2D + H lnt can be written in a 

decoupled form: f(x 2 ) + f(y 2 ), where x = zJi|s| and 
y = zJ±\t\. Since the helical order is related to x and y 
through (axb) oc x 2 — y 2 , we conclude that the collinear 
phase appears if f(x 2 ) has a global minimum at x 2 7^ 0. 
In the quasi-lD case, fix 2 ) — /(0) can be expanded in 
terms of x 2 as Ax 2 + Bx A with 
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1 /A(T) 
A(T) V 2zJ ± 

1 'S- 



A(T) 3 



_ 5+1/2 
86{T) 

2) 128 V 



(1 - 25{T) 2 /if 
1 - 5{T) 2 



where S(T) = A(T)/\(T) (A(T) = A - 2r h cos(q/2)) is 
the renormalized gap. Here we have assumed T 3> A(T). 
Since B is positive for < S(T) < 1, a sufficient condition 
for the collinear phase is A < and a second order phase 
transition to the collinear state occurs at A = 0. Above 
Trvb, S(T) ~ 1 and hence the inequality A < is not 
satisfied for small zJ±. This means Tn < Tbec, where 
Tn is the antiferromagnetic transition temperature. Fur- 
ther lowering the temperature with increasing x, the gap 
collapses to result in BEC of SB. Therefore, we conclude 
Tbec < Tn < Trvb- We have also checked the exis- 
tence of the collinear phase for quasi-2D case by numer- 
ically solving the self-consistent equations without using 
the continuum model. In this way, the instability towards 



the collinear order is a robust feature of the strongly fluc- 
tuating quantum helimagnets, and is essentially different 
from that of classical system with an easy axis anisotropy. 
Now we describe the collinear state a = b = (0, a y ,0) 
(see FigJTJ) . where the 4- fold degeneracy for the energy 
of 7fc<Ta is split into upper and lower branch bands (see 
Figj3]b) . The lower-branch band consisting of linear com- 
binations of 7fc CTa is 2-fold degenerate. The lower branch 
bosons, a k and (3 k , are defined through the Bogoliubov 
transformation as a k = cosh ip k (7^1+ + Cl-ki-)/ V2 — 
sinhipki-jl^ + C^_ kl+ )/V2,/3 k = cosh^ fc (7 fcT _ + 

C*7-fc|+)/^-sinh^ fe (7! i+ + C*7l fcT -)/V2, wh ere C = 
e i7r / 4 and tanh2^ fe = xA(T)/(2fi(fc) 2 - x\{T)) Below, 
we will focus on the low energy dynamics, and neglect 
the upper-branch bosons. This leads to the relation be- 
tween the original bosons b ka : b_ Q/2+k ^ ~ C*&Q/2+fc| 
and b^_ Q/2+kl ~ -(b Q/2+kV 

Helical phase — Next we consider the BEC of the low- 
est modes ao and (3q. This corresponds to the non-zero 
expectation values ofb± Q/2 ^ (\(b Q/2M )\ = \{b_ Q/2l - a )\). 
We obtain the cycloidal helical spin structure as 



S i 

si 

S? 



-sm(Q ■ Ri + n/4)(\(b Qm )\ 2 - \(b Q/2l )\ 2 )/M, 
Scos{Q ■ Ri + 7r/4), 

sin(Q • R t + n/4)((b Qm )(b Q/n ) + c.c.)/Af. (7) 



Here we have used the relaxed constraint J2 ia b\ a bi a = 
2SN '. Now we clarify the relation between the elliptic ra- 
tio and the Bose condensate fraction. If we assume that 
(&q/ 2 |) = while (&Q/2f) 0, Sf becomes zero and the 
elliptic ratio is given by m&/m c ~ \(bQ/ 2 ^)\ 2 /(AfS). In 
this case, the spins are rotating counterclockwise within 
the frc-plane. The clockwise helicity is realized when 
(&q/ 2 |) = while (6q/ 2 |) 7^ 0. Note that the elliptic 
ratio can be much smaller than unity even at zero tem- 
perature due to the strong quantum fluctuation in sharp 
contrast to the classical case. 

Neutron scattering spectra — Now we turn to the neu- 
tron scattering spectra in the helicall phase. For simplic- 
ity, we shall focus on the quasi-lD case with the possi- 
ble relevance to the recent experiment on LiCu202 [H|. 
The magnetic cross section for polarized neutron is given 
by the following correlation functions as (^)±(q) oc 
(SfST q ), where the sign + (— ) corresponds to the paral- 
lel (anti-parallel) neutron spin S n to the a-axis (see Fig. 
[I]). To break the degeneracy of the helicity, we first set 
(^Q/2|) 7^ an d (&Q/21) — 0j i- e -i counterclockwise one. 
We should note here that non-Bragg part is considered 
below, i.e., non-zero q component. In the low energy 
regime, using a k and [3 k , (SQ +q Sl Q _ q ) ~ Fx/q 2 + F 2 /q 
and (S^Q_ q SQ +q ) ~ F 2 /q, respectively, with 



Fi 



\(b. 



Q/2T 



)| 2 



'■A(Tf 
16c 2 
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where q is assumed to be small. The difference 

( S Q+ q S -Q- q ) - ( S -Q-q S Q+ q ) is expressed by the vec- 
tor spin chirality ((Sq +q x S-Q- q ) z ) ji, and is directly 
related to the condensate fraction, i.e., the F\ term. 
The crossover between the F\ and F% terms occurs at 



q c a cx |(&q/2T>| S 



C|| A(T) - u,; ACT)' where a is the 
lattice constant and J|j is the typical energy scale de- 
termined by J\ and J 2 - Another important correlation 
functions, (S q S^: q } (a = x,y), can be observed by the 
S n || c setup. By a similar calculation, one can show 
that {S* Q+q S*_ Q _ q ) = (S*S y _ Q _ q ) for the fluctuating 
part. In the experiment (SQ +q S^Q_ q ) suggests el- 
liptic spiral while (Sq + S"q_ ) indicates circular one. 
This puzzling point would be resolved by our above ana- 
ysis considering the quasi-elasitic component [271 ] . 

Dielectric response — Finally, we examine the dynami- 
cal dielectric response both in the paramagnetic and heli- 
cal phases of the quasi-lD model. Even in the paramag- 
netic and collinear phase, we assume that the fluctuating 
electric polarization is given by Eq. ([T]) [33| . We take the 
mean-field decoupling S, x S 3 = ((&j /1 &j/*)(&jp<r pvhu) - 
h.c.)/(4i) to the ferromagnetic bonds and Si x Sj = 

(( b i fJ , e »P b jp)( b j<? cr a\ e ^ b ™) ~ h-c)/(4i) to the antiferro- 
magnetic bonds. We henceforth focus on the contribution 
from the antifcrromagnetic (J2) bonds since its fluctua- 
tion is stronger than that of the ferromagnetic one. From 
the geometry of the system (see FigJTJ), the polarization 
along the 6-axis P b is always zero. In the paramagnetic 
phase, Im£ aa (o>) =Ime cc (w) due to the rotational sym- 
metry in spin space. The expression for the polarization 
along the a-axis P a is given by 



P a cx (r) 2 /J 2 ) 2^cos(2k x )(ibk<,b- 



ha 



+ h.c). 



For purely ID case, Ime aa {oj) oc (n(uj/2) + 
\/2)/{uj^/uu 2 - 4A(T) 2 ), where n(w) is the Bose distribu- 
tion function. Near the threshold of the absorption, the 
ID van Hove singularity, Ime Q 



M cx l/y/u - 2A(T), 
appears as schematically shown in FigHJa. On the other 
hand, a drastic change of the absorption spectra occurs 
in the helical phase since the low-lying branch bosons 
become gapless (see Figf3]c). In this phase, the energy 
dispersions of the upper and lower branches are given by 
£l + (k x ) — y/c|fc2 -I- 2A(0) 2 and f2_(fc x ) = c\\k x , respec- 
tively. We assume the BEC of SB by the weak interchain 
interaction. The schematic behavior at zero temperature 
is shown in Fig0]b. There are three contributions corre- 
sponding to the processes of two bosons i) in the upper 
branch, ii) in the gapless-lower branch and iii) in both the 
upper and lower branches, respectively. Finally, it should 
be noted that we cannot neglect the one-magnon contri- 
bution coming from the condensed portion in the helical 
phase. This contribution corresponds to that obtained in 
the previous analysis 3J| , but this is much smaller in the 
quantum limit. 




FIG. 4: Schematic plots of Ime a£l (aj) in a) the paramagnetic 
phase and b) the helical phase with T = (A = \/2A(0)). 
Behaviors nearly thresholds are indicated. Blue (solid and 
dotted) lines are the results for purely ID model. Singularities 
are smeared out by the interchain interaction as shown by 
black lines. Insets are schematic boson dispersions in both 
the phases. 
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